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Fractionally charged massive particles (FCHAMPs) appear in extensions of the standard model, 
especially those with superstring constructions. The lightest FCHAMP would be absolutely stable 
and any of them produced during the early evolution of the Universe would be present today. The 
thermal production, annihilation and, survival of an FCHAMP, a lepton L with electroweak (i.e., 
U(1)y) but no strong interactions, of mass ttil and charge Ql (in units of the charge on the electron) 
are explored. The FCHAMP relic abundance is determined by the total annihilation cross section 
which depends on tul, Ql and on the available annihilation channels. Since massive (tul > 1 GeV) 
charged particles (Ql > 0.01) behave like baryons (heavy ions), primordial nucleosynthesis and the 
cosmic background radiation temperature anisotropies limit the FCHAMP relic density. Requiring 
that Ql < Qb /5 leads to a constraint on the Ql - rnL relation. Further constraints on Ql and 
iriL are provided by the invisible width of the Z (Ql > 0.16 for vtil < A/z/2) and by accelerator 
searches for massive, charged particles. Our key result is to exploit the fact that in the early 
Universe, after L freeze-out but prior to e recombination, the negatively charged L~ will combine 
with the more abundant alpha particles and protons to form tightly bound, positively charged 
states (fractionally charged heavy ions). The Coulomb barriers between these positively charged 
L~ a and L~p (L~pp, L~ act, ...) bound states and the free L + suppresses late time FCHAMP 
annihilation in the interstellar medium (ISM) of the Galaxy and on Earth, limiting significantly 
the late-time reduction of the FCHAMP abundance compared to its relic value. The surviving 
FCHAMP abundance on Earth is orders of magnitude higher than the limits from terrestrial searches 
for fractionally charged particles, closing the window on FCHAMPs with Ql > 0.01. However, as 
Ql approaches an integer (e.g., \Ql — n\ < 0.25) these searches become increasingly insensitive, 
leaving some potentially unconstrained "islands" in the Ql - tul plane which may be explored by 
searching for FCHAMPs in the cosmic rays. 



I. INTRODUCTION 

We consider the cosmological and astrophysical vi- 
ability of a fractionally charged massive particle, an 
FCHAMP, the lightest of which is stable on cosmolog- 
ical timescales. While there is an extensive literature on 
integer charged CHAMPs and neutraCHAMPs and their 
neutral counterparts in CHAMP-proton and CHAMP- 
alpha particle bound states (see, e.g., JTH5]), in very few 
of these papers were fractionally charged FCHAMPs ex- 
plored in any detail. Many of these previous constraints 
on CHAMPs do not apply to FCHAMPs. Our study is 
restricted to a fractionally charged lepton L (an L par- 
ticle - antiparticle pair), an 5/7(3) and SU(2) singlet 
with weak hypercharge Y = Ql, which couples to the 
photon and the Z. This choice fixes the total L annihi- 
lation cross section, determining the relic abundance of 
L pairs in the interstellar medium (ISM) of the Galaxy 
and on Earth. Requiring that the relic abundance (mass 
fraction) be bounded limits a combination of the mass, 
ttil, and charge, Ql, providing a lower bound to Q L as a 
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function of m^ |6- 8 . During the post-annihilation freeze 
out evolution of the Universe, prior to recombination, the 
negatively charged L~ will combine with alpha partcles 
and protons to form tightly bound, positively charged, 
heavy ions such as La, Lp, Lpp, Laa. The net posi- 
tive charges of these exotic heavy ions suppresses late 
time annihilation of the bound L~ with the free L + . 
This suppression leads to a terrestrial ratio of FCHAMPs 
to baryons which is so high that terrestrial searches for 
fractionally charged particles [5] challenge their very ex- 
istence. In fjn] the particle physics motivations for the 
existence of FCHAMPs are outlined along with general 
comments on existing experimental limits. In |III| the 
relic FCHAMP abundance is calculated and limits are set 
to rriL and Ql- Constraints on Ql and mt from collider 
searches and from the invisible width of the Z are ana- 
lyzed in §IV| and those from terrestrial searches for frac- 
tionally charged particles in fjVl The post-annihilation 
freeze out capture of negatively charged FCHAMPs by 
protons and alpha particles is described and the ratio 
of Lp and La to baryons in the present Universe is es- 
timated in §VI| In §VII| we compare our results to the 
bounds from the terrestrial searches to decide if they ex- 
clude FCHAMPs - or not. Our conclusions are summa- 
rized in fvTHl 



II. MOTIVATION 

All known color-singlet particles have electric charges 
that are integer multiples of e, the charge of the positron. 
However, it is straightforward to construct theories with 
fractional or even irrational charges. Assuming the stan- 
dard model SU(2) x (7(1) gauge group, the weak hy- 
percharge Y assignments are arbitrarjrl as are the as- 
sociated values of electric charge Q — T| + Y, where 
T 3 is the third generator of weak isospin. Even "or- 
dinary" values of Y can lead to fractional charges for 
the "wrong" SU(3) or SU(2) assignments, such as an 



5(7(3) and 5(7(2)-singlet state with Y 



The hy- 



percharge assignments are no longer arbitrary when the 
standard model is embedded in a simple grand unified 
(GUT) group, such as 5(7(5) [TD], but more complicated 
GUT embeddings allow for fractional charges [IT] . 

Fractionally charged states are extremely common in 
superstring constructions [T2], H3]^] Such fractionally 
charged particles could acquire very large (string-scale) 
masses, or could be confined into integer-charged bound 
states by strongly-coupled hidden sector forces |15rj 
However, they could also occur in the massless sector 
of the theory, e.g., acquiring masses in the TeV range 
from effects related to electroweak physics. 

Another motivation is that kinetic mixing between 
two 17(1) gauge bosons [35] can under some circum- 
stances induce small (e.g., millicharged) and generally 
irrational electric charges for hidden sector particles*] In 
this paper, however, we will only consider larger charges, 
Ql>0.01. 

There have been many searches for fractionally charged 
particles in fixed target and collider experiments, quark 
searches (e.g., Millikan-type experiments), cosmic rays, 
and bulk matter searches, as described in a recent re- 
view [9F] The lightest fractionally charged particle would 
have to be absolutely stable (or at least stable on cosmo- 
logical time scales given the stringent limits on charge 
non-conservation) , and would therefore contribute to the 
cosmological mass density unless it were so heavy that it 
was not produced subsequent to inflation. 



1 Appropriate combinations of non-canonical fermion charge as- 
signments are needed for the cancellation of anomalies. 

2 This is apparent from the Schellekens theorem [14) . which essen- 
tially states that any weakly-coupled heterotic string theory with 
realistic gauge couplings must involve either fractional charges, 
an unbroken SU(5) gauge symmetry, or higher Kafi-Moody em- 
beddings of SU(3) or 5(7(2). Examples with fractional charges 
include both heterotic |15H18| and intersecting D-brane 19 -21 
constructions. For general discussions, see 122112.31 . Recent con- 
structions which avoid fractional charges are described in |24j . 

3 Such bound state cryptons have been postulated as sources of ul- 
tra high energy cosmic rays [25] or as dark matter candidates [261 - 
[28]. 

4 See, e.g, [29H32) . Observational constraints are summarized 
in I33H35] . 

5 Earlier reviews include [36^ 38 . More general reviews of massive 
stable particles include 09 40 . 



Masses smaller than Mz/2 are excluded by the invisi- 
ble Z width unless the charge is small (e.g., Ql < 0.16; 
see i IV ) . There have been many searches for heavy long- 



lived particles at accelerators and colliders, but these are 
difficult to summarize because the analyses are often pre- 
sented for specific assumptions concerning the charges 
and other particle properties. For example, some lim- 
its only apply to fractionally charged particles that carry 
color, while others are presented for specific models, such 
as scalar leptons in supersymmetry. Furthermore, some 
searches are insensitive to particles with charges below a 
critical value such as 2/3 because they are not sufficiently 
ionizing. Searches for fractionally charged particles pro- 
duced in e^- annihilation have excluded charges > 2/3 for 
masses below ~ 102 GeV [SI SI]- There are more strin- 
gent limits from the Tevatron [121411] and LHC [332 S5] 
for heavy quarks or for long-lived integer charged slep- 
tons or charginos. In particular, a DO lower limit of 206 
GeV on the mass of long-lived charged gauginos [33] will 
be translated into a limit on fractionally charged particle 
masses in SJTV 



For definiteness we consider here the simplest exam- 
ple, that of a new, fractionally charged, massive spin- 
1/2 particle^ an FCHAMP "L", that is an 517(3) and 
5(7(2)-singlet with weak hypercharge (and therefore elec- 
tric charge) Y, which therefore couples to both the pho- 
ton and the Z. In calculating the FCHAMP relic abun- 
dance the total ^ annihilation cross section enters. In 
our estimate we consider the dominant annihilation chan- 
nels L + L~ — > //, where / is a standard model fermion 
(including the t for m^ > m t ), as well as the diboson 
channels L+L~ -» 77, 7Z, ZZ, ZH, and WW, where H 
is the Higgs boson (here we assume Mh = 120 GeV). 



III. FCHAMP RELIC ABUNDANCE 

The relic abundance of FCHAMPs may be expressed 
as a ratio of their mass density to the critical mass den- 
sity [7] through the dimensionless mass density parame- 
ter fl L = (pL/pc)o, where pw = (pl++Pl-)o = 2m L n L0 , 
and p c o, the present critical mass (energy) density, de- 
pends on the present value of the Hubble parameter, 
H a = 100ft km Mpc^s- 1 (h « 0.70 [32 HE]), Pco = 
1.05xl0 _5 /i 2 GeV/cm 3 . For the range of L ± annihilation 
cross sections of interest here and, for n l0 = 410.5 cm' 3 , 
corresponding to T = 2.725 K [4"5] , 
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The thermally averaged annihilation rate factor, 
(o"u)ann, in units of cm 3 s _1 , is evaluated using the tech- 



6 The lightest fractionally charged state could also be a scalar, es- 
pecially in supersymmetric theories. The annihilation of charged 
scalars into fermion pairs is P-wave suppressed, implying even 
more stringent cosmological and bulk matter bounds. 



niques in [50] at T — T* « mi/25, the temperature at 
which the abundance of the £ pairs begins to depart 
significantly from their equilibrium abundance 7 1. Note 
that the FCHAMPs are not "frozen-out" at T = T*. For 
T < T* annihilations continue to dominate the evolu- 
tion of L pairs, further reducing their abundance. Only 
when T < T* is the FCHAMP relic abundance "frozen- 
out" , although even then annihilations continue, albeit at 
a rate which is much smaller than the universal expansion 
rate, leaving the relic abundance effectively unchanged. 
The relic ratio of the number of L^ 1 pairs to CMB 
photons is 
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Observations of the CMB temperature anisotropy spec- 
trum constrain the value of the baryon (dominated by 



protons and alphas) to photon ratio, n^/n^ 
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|47j . Barring any significant post freeze-out re- 
newed annihilation (see gV| ) , the present ratio of L^ pairs 
to baryons is 
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In practice, we will evaluate {crv) ann at threshold, i.e., 
T<r t , in part because that is the relevant quantity for 
the subsequent annihilation in the Galaxy and on Earth. 
We have checked that there is little difference (< 13%) in 
the rate factor between T = T» and T <^T* (see Figure 



A. The Annihilation Rate Factor 

The total annihilation rate factor, {av) an n, the ther- 
mally averaged product of the annihilation cross section 
and the relative velocity, may be written in terms of the 
two-photon annihilation rate factor, (<7v) 77 , multiplied 
by N, a function of Ql and rriL (jn L = IOOtoioo GeV), 
which accounts for all the open annihilation channels. 
N r^ (crw) arm /(cru) 77 , where (ctw) 77 = na 2 Q 4 L /m 2 L = 



2.2 x 10~ 25 Q L /m 2 00 cm 3 s _1 (we have adopted a — 
1/128). In computing N we have included, in addition to 
annihilation to two photons (iV 77 = 1), the contributions 
from annihilations to all lepton and quark pairs (with 
m < ttil) along with contributions from annihilations to 
boson pairs WW, WZ, ZZ, ZH, and Z-y, where H is 
the Higgs boson, for m L > (M Bl + M B J/2. (We take 
M H = 120 GeV). N may be written as N = a + b/Q 2 L , 
where a results from annihilation into 77, jZ, and ZZ, 
while b represents annihilations to fermion pairs (three 
generations of quarks and leptons) as well as into WW 
and ZH. With these definitions, 
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1. A Toy Model: Pure QED 

To understand (and anticipate) our results for a = 
a{mL,), b = 6(mi) and f^/1 2 in the standard model (SM), 
it is useful to consider L 1 * 1 annihilation in a "toy model" , 
"pure QED". This model assumes that the FCHAMPs 
interact purely electromagnetically and that there are no 
W or Z bosons (or Higgs). In this model the only anni- 
hilation channels arc 1, — > 77 and li — > / . In this 
case, a,Q ED — 1 and each fermion pair with nif > m^ 
contributes to b® ED an amount oc qj. For three gener- 
ations of quarks and leptons, assuming that m^ > rrib 
(or, m L > m b ), b QED = (20 + 4F(m t /m L ))/3 where, to 
account for the top quark threshhold at mr, = rn t , 
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FIG. 1. Contours of fixed values of Ql in the Ql — vul plane 
for the QED toy model. The red curve corresponds to Q,l = 
Oc (cold dark matter), the blue curve to Ql = Ob (baryons), 
and the magenta curve to Ql = Ob/5; see the text for details. 
The regions above the curves are allowed. 



If n L h 2 is fixed, Eq. 4 leads to a quadratic equation 
for Q L as a function of m 2 ^ which depends on a = a(mi) 
and b = 6(tol). The resulting Q L - m L relation for 
the pure QED toy model is shown for three choices of 
£Il in Figure [l] The small "glitches" at toioo ~ 1-7 
reflect the enhanced annihilation cross section from the 
opening of the top quark channeF] A first, naive choice 



In principle there are also threshold effects for mj, = mi,,m c , m T 



for a lower bound to ft^h 2 is from the requirement that 
the FCHAMP relic density not exceed the density in 
cold dark matter, Q L h 2 < Vl c h 2 w 0.11 07]. The re- 
gion above the red curve in Figure [T] satisfies this con- 
straint. However, charged, massive FCHAMPs do not be- 
have like cold dark matter, they behave more like heavy 
ions or baryons [35) . The region above the blue curve 
in Figure u\ is consistent with fl^h 2 < (^b^ 2 )bbn ~ 
0.022 [51 . However, the relic abundances of the light 
elements, especially deuterium and helium-4, require 
this value of the relic density be provided by baryons, 
so the limit on the FCHAMP density must be even 
smaller (53]. Following the analysis of Dubovsky, Gor- 
bunov and Rubtsov [52], the magenta curve in Figure [l] 
represents fl L h 2 = D, B h 2 /5 « 0.0044. It is the re- 
gion above the magenta curve in Figure [l] that is consis- 
tent with the density constraint; we use this constraint 
(fi^/i 2 < 0.0044, £Il < CIb/S) in our subsequent discus- 
sion. 
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FIG. 2. The contributions to o from annihilation to 77 
(gray), jZ (blue), and ZZ (red), as well as the bosonic con- 
tributions to 6 SM from WW (purple) and ZH (black), as a 
function of ttll = 100 mioo GeV. Note that a^ ED = a 77 = 1. 
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FIG. 3. The contributions to b from annihilation to c, /1, r 
(red), u,c,t (blue), d,s,b (gray), and fi,i = 1,2,3 (purple) 
are shown as a function of niL = 100 mioo GeV. 




etc. These would be almost invisible on the scales of the figures 
and are not explicitly included. Our treatment of final state 
hadrons as free quarks should be a reasonable approximation for 
the ranges (Ql > 0.01, raj, > 1 GeV) considered here. 



FIG. 4. The standard model (solid curve) and QED (dashed 
curve) values for b from L + L~ — >• //, W + W~ , ZH as a func- 
tion of rriL = 100 mioo GeV. 



2. The Standard Model 



With the pure QED toy model as background, we now 
return to the calculations of a, b, and the Q L - m L rela- 
tion for the standard model (SM). In addition to direct 
two photon annihilation, a SM receives contributions from 
annihilation to 7Z and to Z pairs, as shown in Figure [5] 
The individual contributions to 6 from annihilations 
into WW, ZH, and fermion pairs are shown in Figures[2] 
and [3] The total values for 6 SM (solid curve) and & QEt > 
(dashed curve) are shown in Figure |4l and the detailed 
formulae for the cross sections are given in the Appendix. 
Notice the enormous enhancement due to the Z reso- 
nance in the annihilation cross section for mi, ~ Mz/2. 

A comparison between the relic density constraints for 
the toy model (dashed curve) and the standard model 
(solid curve) which follow from the requirement that 
£l L h 2 < 0.0044 are shown in Figure pi The enhanced 
SM annihilation cross section for m L ~ Mz/2 is respon- 
sible for the dip in the Qr, - m L relation at low mr,- 
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FIG. 5. Contours of Ql = ^b/5 in the tul — Ql plane for the 
QED toy model (red dashed curve) and the standard model 
(blue solid curve) using the annihilation rate factors calcu- 
lated at threshold; see the text for details. The regions above 
the curves are allowed. For comparison, the gray dotted curve 
shows the constraint for the QED toy model using the anni- 
hilation rate factor evaluated at T = T* = m_L/25. 



FIG. 6. Constraints in the Ql — rriL plane. The magenta 
curve corresponds to Ql — S7b/5 for the standard model. The 
regions above the magenta curve (Ql < £1b/5) are allowed. 
The region above the red curve is excluded by the constraint 
from the invisible Z-width. The region above the blue line at 
Ql — 2/3 and to the left of the dashed blue line is excluded 
by collider searches. See the text for details. 



IV. ACCELERATOR AND Z- WIDTH 
CONSTRAINTS ON Q L AND m L 

The review article by Perl, Lee and Loomba [9] pro- 
vides a summary of the searches for fractionally charged 
elementary particles. The best accelerator/collider con- 
straint for a generic color-singlet fermion comes from the 
OPAL collaboration gT] at LEP 2, excluding Q L > 2/3 
for m L < 102 GeV at 95% c.l. 

The DO collaboration at the Tevatron has set 95% c.l. 
lower limits of 206 GeV and 171 GeV, respectively, on 
the masses of long-lived charged gauginos and Higgsi- 
nos [23] • For m L < 300 GeV these can be approxi- 
mately translated into limits on unit-charged FCHAMPs 
by combining the theoretical NLO cross sections and ex- 
perimental limits (as a function of mass) presented by 
DO with the parton-level FCHAMP and gaugino Drell- 
Yan cross sections given in the Appendix. We find that 
m L < 142 GeV for Q L — 1. The result can be extended 
to other values of Q L if one makes the additional assump- 
tion that the acceptance depends only weakly on Q L . We 
assume that this is valid for Q L > 2/3, corresponding to 
m L < 115 GeV, in which case the combination of the 
OPAL and DO collider limits excludes the region in the 
upper left corner of Figure [6] bounded by the solid and 



dashed blue lines. 

The Z may decay to an FCHAMP pair provided that 
iUl < Mz/2. The limit to the invisible width of the Z, 
often expressed as a limit to the number of equivalent 
neutrinos, N v = 2.984±0.009 [5"3" .provides a lower limit 
to Ql(> 0.16) as a function of mi ®| This bound is shown 
as the red curve in Figure [6] 

The regions above the magenta curves, below and to 
the right of the red curve, and below and to the right of 
the blue line and dashed curve in Figure [6] identify the 
regions in the Ql - tul plane that are consistent with the 
relic density, collider, and invisible Z-width constraints. 
We next consider the constraints imposed by terrestrial 
searches for fractionally charged, massive particles. 



V. CONSTRAINTS FROM TERRESTRIAL 
SEARCHES FOR FCHAMPS 

The null results of searches for fractionally charged 
particles in bulk matter on Earth (E) [54J, [55] or me- 
teoritic material (e.g., |54[ I56j ) are summarized in the 
review by Perl et al. |5]. Their 95% confidence upper 
limits for terrestrial material from the oil drop or mag- 
netic levitometer techniques, (n L /n B ) E < 1.3 x 10 -23 , 
and for meteoritic material (M), (n^/nr^M < 7. Ox 10 -22 , 
provide enormously stringent limits to free terrestrial 
FCHAMPsf] Since the FCHAMP density is p L = 
2m^rii and the baryon density is p B — ^b^B; the ratio 
of relic FCHAMPs to baryons by number may be written 
as 



and, 



n L _ 9.4 x 10" 
n B mioo 



Q, 



4.7 x 10" 3 



Tiioo 



n B /5j Ql(aQ 2 L + b) 



(6) 



The relic abundance of FCHAMPs is many orders of mag- 
nitude larger than the upper bound to their density set 
by the terrestrial constraints. However, long after anni- 
hilation freeze out in the early Universe, when the solar 
system and Earth form, annihilation resumes in the re- 
sulting high density, low temperature environment. In- 
deed, at the low terrestrial temperature (T ~ 300 K), the 
annihilation cross section is "Sommerfeld enhanced" [57] 
from its early Universe value, 



(7) 

7.2 x 
10 4 Q 2 L m 10Q , (<71>)e is evaluated at threshhold (T B <C T*) 
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s(v) 



'I 



1 - e~i i7»l 



>T). 



(8) 



For equal numbers of positively and negatively charged 
FCHAMPs, renewed annihilation on Earth is rapid pro- 
vided that T® nn At = n L (av) E At > 1, where At « 
4.5 Gyr and n L = n B (n L /n B ) « 6.4 x 10 23 (ti l /?ib) [H]- 
In this case, the annihilation on Earth efficiently reduces 
the relic L/B ratio to a present value of 

6.9 x 10~ 22 m^ 



n L 



1 



nB/ E n B {<Jv) E At Q L {a 



b) 



The requirement that (nL/n B ) B < 1.3 x 10 23 leads to 
the constraint on Ql as a function of raj, , 

Qi(aQl + b)>53m%l (10) 

which is shown by the dashed, black curves in Figure [7J 
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FIG. 7. As for Figure [6] constraints in the Ql — ttil plane. 
The black, dashed curve corresponds to an allowed terrestrial 
FCHAMP abundance resulting from L^ annihilation over the 
Earth's history, provided that all relic L~ were free (i.e., 
not bound to alphas or protons). The allowed region is above 
the black, dashed curve (subject to the collider (blue) and 
invisible Z- width (red) constraints). See the text for details. 



This is the 95% c.l. limit, taking into account that the physical 
value of N„ — 3 is positive semi-definite. The total Z- width 
obtained directly from the lineshape gives a somewhat weaker 
lower limit, Q L > 0.24 for m L < M z /2. 

On average, the net FCHAMP charge on a drop would be zero, 
but there are fluctuations in the numbers of L+ and L~ in any 
drop (including those bound to protons and alphas) and many of 
the fluctuations should yield drops whose net charge differs from 
an integer. 



The situation is quite different for free FCHAMPs in 
the ISM where n B ~ 1 cm" 3 and T sa 10 4 K. In this 
case, 7/ism w 1-25 x 10 4 Q 2 (m^/T^ 1 / 2 (T = 10 4 T 4 K) 
and the ISM annihilation rate for free FCHAMPs, 



^ISM 



= n L (av)isu~ 1.3x10 23 
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noo 



«B 



T 



1/2 



(11) 



is far too small to deplete any free L ± pairs in the Galaxy 
over the age of the Universe (t — 14 Gyr, T l ®j£t < 1). In 
contrast to free FCHAMPs on Earth, those in the ISM 
are effectively immortal. 



A. Constraints on Ql And mt 

As may be seen from Figure [7J for Ql < 2 and 
m,L < 300 GeV, the only potentially allowed values of 
Ql and tul are found above the dashed black curves, 
subject to the collider (blue) and invisible Z-width (red) 
constraints. Although the regions above the relic abun- 
dance curve (magenta) are consistent with VIl < ^b/5, 
in the regions below the dashed black curves the terres- 
trial abundance of FCHAMPs, (n£/nB)E, is too large. 
For example, over the entire range of Ql and rriL shown 
in Figure [71 in the absence of further annihilation, the 
relic FCHAMP to baryon ratio exceeds ~5x 10~ 7 . In 
the regions above the dashed black curves annihilation 
on Earth can reduce this ratio to < 1.3 x 10~ 23 , consis- 
tent with the terrestrial constraints. However, if some 
FCHAMPs, even a tiny fraction of them, are insulated 
from annihilation, even these allowed values of Ql and 
TCIl values may be excluded by the terrestrial constraints. 
We address this issue next. 




FIG. 8. The fraction of negatively charged FCHAMPs, f L , 
captured into bound states with alphas (red), with protons 
(solid blue) and with two protons (dashed blue) as a function 
of Ql ■ The magenta curve shows the sum. 



VI. BOUND STATES OF IT WITH ALPHA 
PARTICLES AND PROTONS 

As the early Universe expands and cools, annihilation 
of free Li^ pairs effectively ceases but, at sufficiently low 
temperatures in the early Universe, after primordial nu- 
cleosynthesis (big bang nucleosynthesis: BBN) has con- 
verted some nucleons into alpha particles, collisions be- 
tween alphas and protons and the negatively charged 
FCHAMPs (L~) produce bound L~a, and L~p systems 
(= Lp, La) [2HS] with charges 2 — Ql and 1 — Ql respec- 
tively For Ql > 1 (2), Lpp (Lena) bound states with 
electric charges 2 — Ql (4 — Ql) will form, replacing the 
Lp (La) bound states. Because of the exceedingly large 
ratio of cosmic background photons to FCHAMPs, these 
bound states can't form until T <C £*, typically not un- 
til T < £"h/40, where E^ is the binding energy of the 
system (or, of the system in its first excited state) [55] , 
Since for the La and Lp bound states this occurs after 
BBN has ended (Tbbn ^ 30 keV), these bound states 
will not affect the BBN-predicted relic abundances (un- 
less Ql > 2). However, since there are considerably 
more baryons (nucleons) than FCHAMPs, a significant 
fraction of all L~ will emerge from the early Universe 
"dressed" by alphas or protons in tightly bound, posi- 
tively charged, heavy ions. Collisions between the free 
L + and these positively charged heavy ions are "Gamow 
suppressed" |59j by their mutual Coulomb repulsion (the 
flip side of Sommerfeld enhancement between oppositely 
charged particles). 



A. Capture of L by Alpha Particles 

After BBN has ended the fraction of all baryons in 
alpha particles is n a /n B = Y P /4 m 1/16 [5T], where 
Yp is the primordial helium mass fraction. Using the 
estimate in eq. 6 for the baryon to FCHAMP ratio, 



n L _ 0.015 / n L 
n a mioo \^b/5 



< 



0.015 

TOioo 



(12) 



In the early Universe alpha particles outnumber 
FCHAMPs by a large factor. At sufficiently low temper- 
ature, T < Et,(La)/40, free L~~ will combine with alpha 
particles to form hydrogen-like bound states. For Ql < 
2, the La bound state has a net positive electric charge 
2 - Ql > 0, a binding energy E b (La) s=a 400Q| keV, and 
a ground state radius r(La) « 3.6/Ql frd_] This bound 
state forms when T < 10Q| keV (< 40 keV for Q L < 2), 
after the alpha particles have formed during BBN. 

However, just because there are more alpha particles 
than FCHAMPs and La bound states can form, does not 



10 Since the radius of the alpha particle is comparable to the ground 
state radius of the La ion, the L~ will not see the full charge of 
the alpha particle and there are small corrections to the binding 
energy and radius which are ignored here. The difference between 
the alpha particle mass and the reduced mass of the La system 
has also been ignored. 



guarantee that they will form. An estimate of the frac- 
tion of negatively charged FCHAMPs captured by alpha 
particles may be found by comparing La formation to 
ordinary proton - electron (re)combination, which oc- 
curs when the temperature drops below ~ Eb(H)/AQ « 
0.3 eV. At this time the recombination rate, T pe — 
n e {av) pe , greatly exceeds the universal expansion rate, 
H, (T pe /H) rec « 2.6 x 10 3 [55]. The La (re) combination 
rate, ^La.rec = (n a (v<j) L a ) rec may be scaled to the pe 
recombination rate, since 



1 
12 



T 



2.1 x 10 12 Q|, (13) 



where T La ,reclT pe ^ ec = E b (La) / E b (pe) 
4Q 2 L (m a /m e ) « 2.9 x 10 4 Q|, and 



\M)Lc 



{(TV 



pe 



(2Ql) 



3/2 , rp \ 1/2 



= C2Ql) ( 



T, 



pe 



so that 
'Tl. 

r 



pe 



(2Ql) 6 (T l . 
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T, 



pe 



Tlci 



7.4 x 10~ 8 Q|, (14) 



1.6xl0 5 Q|. (15) 



Since pe (re)combination occurs close to when the Uni- 
verse is making the transition from radiation dominated 
to matter dominated, 

(H pe /H ) rec = VLn(l+ Z pe ^ rec ) +f2M(l+%>e,rec) • (16) 



^pe,rec 
pe,rec 



For Qm/^r = 1 + z eq « 3.2 x 10 3 g7] and 1 
1.1X 10 3 , H pe , rec ~ 4.9x 10~ 14 s" 1 , so that (T/H)pe 
2.6 x 10 3 . In constrast, La (re)combination occurs when 
the Universe is completely radiation dominated so that 

H La ,rec/H = il^il+ZLa.rec) 2 , where (1 + ZLa.rec) / 0- + 
Zpe,rec) — -*-Lot,rec 



/J-pe,rec: SO that 



(T La /H) 



10.6Q| 



J pe 
Tha 



and {T La /H) rei 
» 1, (T La /H) r 



k 3.6 x l(T 4 Qt, 

(17) 
i 0.94Q£. Unlike (T/H) pe , rec which is 
increases from > 9.4 x 10~ 9 to < 15 
as Ql increases from > 0.01 to < 2. An estimate of 
the fraction of FCHAMPs captured by alpha particles is 
given by f La = [1 - exp(-(r La /_ff) rec )] [3] which, for 
0.01 < Ql < 2, increases by some eight orders of mag- 
nitude from ~ 10 -8 to 1.0; see Figure [8j Although for 
parts of this Ql range most of the negatively charged 
FCHAMPs may escape capture by alpha particles leaving 
them vulnerable to annihilation, a non-negligible fraction 
of them will be sequestered in tightly bound, positively 
charged La ions, insulating them from annihilation with 
their free, positively charged counterparts. This leads to 
a ratio of the surviving FCHAMPs to baryons on Earth, 
/l(^l/^b)j far in excess of the limits set by the terres- 
trial searches (< 1.3 x 10 -23 ) [§]. 



While the evolutionary history of meteoritic mate- 
rial is uncertain, such samples may represent a bridge 
from interstellar to terrestrial matter. Since in the ISM 
r ann At < 1, even for free FCHAMPs, the interstellar 
FCHAMP to baryon ratio is the relic value (see eq. 6), 
far exceeding the limits set by searches for fraction- 
ally charged particles terrestrially or in the meteorites 
(< 7.0 x 10- 22 ) 0. 



B. Capture of L By Protons 



Since the fraction of all baryons in the early, post-BBN 
Universe which are free protons is n p /n^ = 1— Yp « 
0.75 [5J, the FCHAMP to proton ratio, by number, is 



n L 



1.25 x 10~ 3 / fl L \ 1.25 x 10~ 3 



mioo V^b/5 



< 



m W o 



(18) 



The Lp bound state forms later than the La ion, when 
T < T L p,rec ~ 0.5Q| keV. Any L~ which had failed 
to form La ions may combine with protons to form Lp 
bound states at this lower temperature. For Ql < 1, the 
Lp ion has a net positive electric charge 1 — Ql > 0, a 
binding energy Eb(Lp) ss 25Q| keV, and a ground state 
radius r(Lp) ~ 29/Ql fin- 
As was the case for the capture of L~ by alpha par- 
ticles, it must be asked what fraction of any free L~ , 
those which escaped earlier capture by alpha particles, 
will be captured by protons when T < TL P . re c- Scaling 
the above analysis for La to the Lp system results in 
(TL P /H) rec m 2A3Q L . The corresponding Lp capture 
fraction, 

f L p{i - f La ) = (i - cxp(-(r/i7) Lp . r . ec ))ex P (-(r/ J ff) LQ ) 

(19) 
is shown in Figure^ (solid blue curve). 

The fraction of all negatively charged FCHAMPs 
which are either captured by protons or by alpha par- 
ticles, 

fhtot = 1 - exp(-(r/ J ff) Lp!rec - (T/H) La ^ ec ) 

= l-cxp(-3.37Qi), (20) 

ranges from > 5.4 x 10~ 3 to < 0.97 as Ql increases from 
0.2 to 1 (see Fig. [8]). 

It is likely that later in the evolution of the Universe, 
e.g., in the galactic interstellar medium or, perhaps in 
stars, strongly exothermic charge transfer reactions Lp + 
a — > La + p will rearrange all the bound, negatively 
charged FCHAMPs into La heavy ions with charge 2 — 
Ql>1- 



/. Lpp Bound States 

For Ql > 1, the Lp bound state is negatively charged, 
which would result in an enhancement of the annihila- 
tion with free L + . However, when Ql > 1, two-proton 
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bound states with L~ become possible. Scaling to the 
binding energy of the H~ system (replacing the proton 
with an L~ and the two electrons with two protons), for 
Ql > 1, Eb(Lpp) > 1.4 keV. Repeating the above anal- 
ysis for the formation of the Lpp bound state leads to 
the dashed blue curve in Figure [HI For Ql > 1, the Lp 
bound states are "replaced" by Lpp bound states with 
positive charge 2 - Ql (for Ql < 2). 



2. Laa Bound States 

The analysis here has been limited to Ql < 2. For 
Ql > 2 the La bound state has negative charge and this 
might lead to an enhancement in the late time annihila- 
tion of L pairs. However, in analogy with the formation 
of the Lpp bound state when Ql > 1, a two-alpha bound 
state can form with L~ when Ql > 2. Scaling to the 
binding energy of the neutral helium atom (replacing the 
helium nucleus with an L" and the two electrons with 
two alpha particles), for Ql > 2, Eb(Laa) > 180 keV. 
For Ql > 2, the La bound states would be "replaced" 
by Laa bound states with positive charge 4 - Ql (for 
Ql < 4). We note that for Ql > 2, La (re) combination 
would occur for T rec > 40 keV, during BBN (after the 
formation of alpha particles). It is not unlikely that 
these negatively charged particles might alter the BBN 
abundances of nuclides heavier than helium by helping to 
bridge the gaps at mass-5 and mass-8, perhaps leading 
to the BBN formation of L~ bound states with C, N, O 
nuclei which would have even larger positive charge. In- 
vestigation of this possibility is beyond the scope of this 
paper. 



VII. DISCUSSION 

Consider the consequences of the capture of L~ by al- 
pha particles and protons in the early Universe. Over 
the range 0.01 < Q L < 2/3 (1 < m L < 115 GeV), 
the fraction of the relic, negatively charged FCHAMPs 
in positively charged bound states with alpha particles 
or protons increases from ~ 3 x 10~ 8 to ~ 0.49. Even 
though this leaves most of the relic L~ free to annihilate 
on Earth with the relic L + , the bound L~ are insulated 
from annihilation. For vul < 115 GeV, the surviving 
ratio of FCHAMPs to baryons on Earth, /L(ni/«B), ex- 
ceeds ~ 9 x 10 -11 , more than 12 orders of magnitude 
higher than the terrestrial upper limit. All combinations 
of Ql and tul in this regime are excluded. The conflict is 
even worse at higher masses. For 142 < niL < 300 GeV, 
the DO constraint jH] provides an upper bound to the 
charge which increases with niL, from QL,max — 1-0 at 
m L = 142 GeV to Q L ,max = 10.3 at m L = 300 GeV, 
leading to a lower bound to the relic FCHAMP to 
baryon ratio which decreases from n^/riB > 6 x 10~ 4 
at m L = 142 GeV to n L /n B > 7 x 10~ 7 at m L = 300 
GeV. For tul > 300 GeV there is no upper bound to 



Ql so that for fixed tul the lower bound to the relic 
FCHAMP to baryon ratio decreases as a power of Ql, 
{nL/n-B)min oc Q^ 4 - However, to reduce the FCHAMP 
abundance on Earth to a level consistent with the ter- 
restrial upper bound would require a value of Ql far in 
excess of the range considered here. 



A. Prying Open The Window For FCHAMPs 

The very strong conclusions reached here, precluding 
the existence of FCHAMPs with Q L > 0.01 (m L > 
1 GeV), follow from the very low upper bounds to the 
FCHAMP to baryon ratio on Earth inferred from many 
oil drop and magnetic levitometer searches using ter- 
restrial material (see [9] and references therein). How- 
ever, these searches become insensitive to FCHAMPs 
as Ql approaches an integer. According to Perl et 
al., (Oj, FCHAMPs may not have been excluded if 
\Ql — n\ < 0.25, where n is an integer. Islands in the 
Ql ~ m L plane satisfying \Ql — n\ < 0.25 and consistent 
with the relic density (fiz, < f2s/5), accelerator, and Z 
width constraints are shown in Figure [9] 




FIG. 9. Islands in the Ql - mL plane consistent with the 
relic density, accelerator, and Z width constraints, for which 
\Q L -n\ <0.25 (n = 0,l,2). 



If, indeed, the combinations of Ql and vul shown in 
Figure [9] have not been excluded by searches in terres- 
trial material, such FCHAMPs may be observable in the 
cosmic rays. The cosmic rays result from the acceler- 
ation of material in the ISM of the Galaxy and if the 
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interstellar FCHAMPs are efficiently acceleratecpj the 
resulting cosmic ray FCHAMP to nucleon ratio might be 
detectable. Since free, as well as bound, ISM FCHAMPs 
avoid annihilation, the predicted distribution of fraction- 
ally charged particles in the cosmic rays provides a smok- 
ing gun for FCHAMPs. For example, depending on the 
value of Ql, there will be free L ± with q = ±Ql, along 
with La with q = +(2 — Ql) an d Lp with q = +(1 — Ql) 
(and/or Lpp with q = +(2 — Ql) for Ql > 1 and Laa 
with q = +(4 — Ql) for Ql > 2). The relative abun- 
dances of these charge states is determined by the rel- 
ative capture fractions, f La , f Lp (hpp, /tea), and f to t 
(see Figure [8]). 

To illustrate the potential of the cosmic rays to pro- 
vide a smoking gun for FCHAMPs, consider the following 
example. Suppose that Ql = 5/6 (\Ql — 1| = 1/6 < 
0.25) and tul = 110 GeV. For this choice of param- 
eters, n L = 0.17tt B < Q B /5 and the relic FCHAMP 
to baryon ratio by number is n^/ne = 7 x 10~ 4 . For 
Ql = 5/6, the fraction of negatively charged FCHAMPs 
that are free (q = —5/6) is 1 — f to t = 0.20 so that 
iV(-5/6)/iV(+5/6) = 0.20, N(+7/6)/N(+5/6) = f La = 
0.36 and, N{+l/6)/N{+5/6) = f Lp = 0.44. 



VIII. CONCLUSIONS 

FCHAMPs with Q L > 0.01 and m L > few GeV be- 
have like heavy baryons [33], leading to an interesting 
constraint on their relic mass density [52], Ql ^ ^b/5 
(fi^/i 2 < 0.0044). However, if this constraint were satu- 
rated the relic FCHAMP abundance on Earth would be 
orders of magnitude larger than the limits set by the 
negative results of searches for terrestrial, fractionally 
charged particles. But, if the relic terrestrial FCHAMPs 
consisted of "free" L ± pairs, renewed, Sommerfeld en- 
hanced, annihilation would have occurred over the past 
4.5 Gyr in the high density, low temperature environment 
of the Earth, reducing the relic abundance and resulting 
in a range of Ql - mL values {Ql > 0.2, rriL > Mz/2; 
see Figure UK consistent with the upper bounds inferred 
from the terrestrial searches. However, as we have noted 
(see ( VI I, in the early Universe after L ± freeze out, nega- 



tively charged FCHAMPs will be captured into positively 
charged bound states with alpha particles and protons, 
suppressing any further, late time annihilation. This 
would appear to close the window on any FCHAMPs 
with 0.01 < Q L < 2 (< 4). Although a "natural" upper 
bound for the range of the electric charge of an FCHAMP 
might be Ql < 1, we have extended the analysis here to 
Ql < 2 (or, to Ql < 4 by including Laa). We expect 
that larger values of Ql would also be excluded by the 



11 Note that FCHAMPs have a much smaller charge-to-mass ratio 
than do the ordinary nucleonic cosmic rays which may affect the 
acceleration efficiency. 



negative results of searches for fractionally-charged par- 
ticles on Earth because of bound states with multiple a's 
and/or heavier nuclei. In that case, Eq. [6] would suggest 
that only enormous values of Ql(> 10 4 — 10 5 ) would be 
allowed. We have not considered the very large Ql sce- 
nario here because it would likely require modifying BBN 
and, furthermore, perturbative calculations would break 
down for aQl/4-rr = (9(1). Note that for m L < 300 GeV, 
the DO constraint [12] suggests that QL,max ;$ 10.3, re- 
sulting in a lower bound to the relic FCHAMP to baryon 
ratio, n^/riB >7x 10~ 7 . 

Having argued that the extremely low upper bound 
from terrestrial searches for FCHAMPs has closed the 
window on such particles, we then noted that such 
searches may have been insensitive to FCHAMPs whose 
charge is close to an integer (\Ql~ n\ < 0.25) and we iden- 
tified some islands in the Ql - rriL plane (see Figure [9| 
where fractionally charged massive particles may not 
have been excluded. FCHAMPs with these charges and 
masses would reside in the ISM of the Galaxy and may 
have been incorporated into the cosmic rays. Searches for 
such particles and, if there is evidence for them, the ratios 
of different charge states (±Ql, +(2 — Ql), +(1 — Ql), 
etc.) can provide a smoking gun for the existence of 
FCHAMPs. It is also anticipated that the collider lim- 
its, especially for large Ql and tul, will be considerably 
strengthened at the LHC. 
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Appendix: Annihilation Cross Sections: Standard 
Model 



Here we consider annihilation into // and W + W~ by 
s-channel 7 and Z exchange, into ZH by s-channel Z 
exchange, and into 77, jZ, and ZZ by t- and it-channel 
L exchange. The results below and near the Z pole use 
the Breit-Wigner propagator (s — M§ + iAIzTz) 
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1. Annihilation to // 



3. Annihilation to ZH, where H is the Higgs boson 



Define yf L (yf R ) as the weak hypercharge of the left 
(right) chiral projection of fermion /, so that qf — 
1 )l + V/l = Vf R (y/z, = V 6 for the quarks and -1/2 
for the leptons). The U(l)y boson couples with strength 
g'yf L to /l, with strength g'yf R to fa, and with strength 
g'yh — g'Qh to L. The full expression, including rrif, 
Mz, and Tz is too complicated to give here, but is in- 
cluded in the plots. In the limit m^ 3> to/, Mz, 



{° v )ff 



C f g A Qitm^e w {y} L +q 2 f 
32ttto| 



(A.l) 



where g' = g tan 9\\? and C/ = 3 (quarks) or 1 (leptons) 
is the color factor. Note that e = gsin9w- Eq. |A.1| is 
equivalent to the QED expression 



, nQED C f Tra 2 Q 2 L q 2 f 
M// = —2 • 



(A.2) 



except that (in this limit) it is the U(1)y boson B that is 
actually exchanged. That is, for tol much larger than the 
electroweak scale, one can think of W 3 and B exchange 
rather than 7 and Z exchange. However, the W 3 does 
not couple to L since it is an SU(2) singlet. Including a 
final mass to/ < mi, 



C f g i P f Q 2 L taii i e 



w 



x (6q f r 2 f y fL + (4 - r 2 ) y% + qj (4 - r 2 )), (A.3) 



128ttto 2 



where r/ = rrif/m^ and /3/ = \/l — r*, generalizing 

Eq.JSJ 

The corresponding values of 6?¥ for three families are 
shown as a function of mi in Figure [31 including the 
effects of Mz and T^. 



2. Annihilation to 77, 7Z, and ZZ 



The cross sections for these channels are 
/ s 9 4 Qi sin4 e w 

( CTW )77 = 1( . 2 



(ctv) 



(<™) 



7 Z 



ZZ 



1677771^ 

ff 4 Qi (4 - r|) sin 4 % tan 2 0^0(1 - r z /2) 
ff 4 Qi (1 - r%) 3 ' 2 sin 4 6> w tan 4 9 W 9(1 - r z ) 



4ttto 2 (2 - r§) 2 



(A.4) 



where rz = Mz/rriL. For tol 3> M^, the sum reduces 
to 

.g 4 g|tan 4 6»H/ ,. , 

(g ^ = 16.M 2 ' (A - 5) 

which is just the BB cross section since g 4 tan 4 9w = <7 ■ 



The cross section for the ZiJ annihilation channel is 



{av) 



ZH 



_ g^kfQl tan 4 9 W 

25()Trm 3 L 
16 + 40r| - 2r| (4 + r|) + r^ 



16-8r|+r 4 +r| 7 2 
where r# = M H /m L , 72 = T z /m L , and 

^4to 2 - (A/ z - Mh) 2 ^4to 2 - (M ff + M z ) 



(A.6) 



k f = 



4m L 



For tol ^$> Mz , Mff this reduces to 



{av) 



ZH 



g 4 Q 2 L tan 4 9w 
2567tto| 



(A.7) 
(A.8) 



Annihilation to W + W 



The cross section is 



{av) 



WW 



g^QiPfSin^w 



(4 + 16? 



w 



QAinn 2 L 
3r 6 



17^) 



z (4 + 7§) 



r^(16-8r|+r| + r| 7 l) 



(A.9) 



where r^ = Mw/rriL, rz = Mz/rriL, and /3/ = 
v/l-r^. For to l > M w , M z 



(av)ww 



g A Q\ tan 4 #vk 
2567rm? 



(A.10) 



It is at first surprising that this is nonzero. For large 
TOi the 7 and Z contributions cancel, leading to a factor 
M§/to 2 — > in the amplitude (i.e., one expects only the 
linear combination corresponding to the B to survive, 
but the B doesn't couple to W + W~). However, the po- 
larization vectors for the longitudinal W ± each vary as 
ttil/Mw in that limit, leading to a finite amplitude. The 
same result can be derived using the equivalence theo- 
rem [60j at high m^. 



Annihilation total for two bosons 



The contributions of the diboson channels to 



,SM 



and 

b are shown as a function of TOi in Figure |2j including 
the effects of Mz and Tz- The Higgs mass is taken to be 
120 GeV. 



6. Drell-Yan production of FCHAMPs and 
gauginos 

The high-energy parton-level cross section for 
FCHAMP (L+L ) and charged gaugino (W + W~) 
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production are given by 

<? 4 Qf tan 4 w (yl +<?!)£/ (3-£/) 



&LL = 



'WW 



1447TS 



(<&-yiJ a /»/(3-$) 

1447TS 



(A.ll) 



respectively, where g^ and y^ are the electric charge and 
(left-chiral) weak hypercharge of the initial parton, and 
P f = y/1 - Am 2 js is the final velocity of the FCHAMP 



or gaugino (of mass to). In Eq. A.ll we have neglected 
the parton mass and assumed s S> M%, in which case 
the production proccdcs by s-channel B and W 3 ex- 
change, respectively. The ratio <Jll/<?ww IS therefore 
found to be - 0.17Q| (0.05Q|) for uu (dd). We have 
verified numerically that these ratios continue to hold to 
a good approximation when the full expressions involving 
Z and 7 exchange are utilized. Since uu dominates at the 
Tevatron, we assume that the ratio of pp cross sections 
for FCHAMP and gaugino production is 0.15 Q\ in the 
Ql — m L range relevant to the limits described in jjIV| 
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